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J  Dynamic  fracture  processes  in  structural  materials  are  often  described  in  terms  of  the 
relationship  between  a  measure  of  the  crack  driving  force  and  the  crack  tip  speed.  In  this 
article,  ongoing  research  directed  toward  establishing  a  basis  for  such  a  relationship  in  terms 
of  crack  tip  plastic  fields  is  described.  In  particular,  the  role  of  material  inertia  on  a  small  scale 
as  it  influences  the  perceived  fracture  resistance  of  a  rate  independent  material  is  discussed. 
Also,  the  influence  of  material  strain  rate  sensitivity  on  the  development  of  crack  tip  plastic 
deformations,  and  the  implications  for  cleavage  propagation  and  arrest  in  a  material  that 
can  undergo  a  fracture  mechanism  transition,  is  considered.  The  discussion  is  concluded 
with  mention  of  a  few  outstanding  problems  in  the  study  of  dynamic  fracture,  including 
some  recent  experimental  evidence  that  the  traditional  crack  tip  characterization  viewpoint 
of  fracture  mechanics  may  be  inadequate  under  very  high  rate  loading  conditions. 


KEYWORDS:  dynamic  fracture,  dynamic  crack  propagation,  elastic-plastic  fracture,  high 
strain  rate  fracture,  crack  arrest 
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1.  INTRODUCTION 


Consider  growth  of  a  crack  in  an  elastic-plastic  material  under  conditions  that  are  es¬ 
sentially  two  dimensional.  The  process  depends  on  the  configuration  of  the  body  in  which 
the  crack  grows  and  on  the  details  of  the  applied  loading,  in  general.  However,  if  the  region 
of  active  plastic  flow  is  confined  to  the  crack  tip  region,  and  if  the  elastic  fields  surround¬ 
ing  the  active  plastic  zone  are  adequately  described  in  terms  of  an  elastic  stress  intensity 
factor,  then  it  is  commonly  assumed  that  the  prevailing  stress  intensity  factor  controls  the 
crack  tip  inelastic  process.  With  this  point  of  view,  the  stress  intensity  factor  provides  a 
one-parameter  representation  of  the  input  into  the  crack  tip  zone.  The  viewpoint  mimics  the 
small-scale-yielding  hypothesis  of  elastic-plastic  fracture  mechanics  [1]  but  the  basis  for  it  in 
the  study  of  rapid  crack  growth  is  less  well  established. 

Experimental  data  on  rapid  crack  growth  in  elastic-plastic  materials  is  commonly  inter¬ 
preted  on  the  basis  of  an  extension  of  the  Irwin  crack  growth  criterion  in  those  cases  in  which 
a  stress  intensity  factor  fields  exists.  If  the  applied  stress  intensity  factor  is  Ka  then  a  com¬ 
mon  constitutive  assumption  is  that  there  exists  a  material  parameter  or  material  function, 
say  depending  on  crack  speed,  and  possibly  on  temperature,  such  that  the  crack 

grows  with  Ka  =  Km.  Indeed,  in  the  jargon  of  fracture  dynamics,  such  a  condition  provides 
an  equation  of  motion  for  the  position  of  the  crack  tip  as  a  function  of  time. 

Rapid  crack  growth  in  metals  subjected  to  quasistatic  loading  or  stress  wave  loading  of 
modest  intensity  seems  to  follow  this  constitutive  assumption  to  a  sufficient  degree  so  that 
a  systematic  study  of  some  of  its  physical  underpinnings  is  warranted.  Thus,  in  recent  years 
considerable  effort  has  been  devoted  to  developing  models  to  explain  the  reasons  why  I\'„, 
depends  on  crack  speed,  and  possibly  on  temperature,  as  it  does  for  real  materials.  The 
approach  is  quite  straightforward.  It  is  assumed  that  a  crack  grows  at  some  speed  v  in 


an  clastic-plastic  or  elastic- viscoplastic  material  under  the  action  of  the  input  A„  applied 
remotely  from  the  crack  tip  region.  The  potentially  large  stresses  near  the  edge  of  the 
crack  are  relieved  through  inelastic  deformation  in  an  active  plastic  zone,  and  a  permanently 
deformed  layer  is  left  in  the  wake  of  the  active  plastic  zone  along  each  crack  face  as  the  crack 
advances  through  the  material.  A  solution  of  this  problem  is  then  obtained  for  arbitrary  I\.a 
and  v  in  the  form  of  stress  and  deformation  fields  that  satisfy  the  field  equations  in  some 
sense.  With  this  solution  in  hand,  a  crack  growth  criterion  motivated  by  the  physics  of 
the  process  may  be  imposed  on  the  solution  to  yield  a  relationship  between  Ka  and  v  that 
must  be  satisfied  for  the  crack  to  steadily  advance.  This  relationship  is.  in  fact,  the  material 
function  I\m(v,T)  for  the  model  problem. 

The  way  in  which  I\m  depends  on  t\  T  and  other  system  parameters  depends  critically 
on  the  details  of  the  fracture  separation  process,  that  is,  whether  it  is  a  void  nucleation 
and  ductile  hole  growth  mechanism  or  a  cleavage  mechanism,  whether  there  is  a  strain  rate 
induced  elevation  of  flow  stress  or  not,  whether  the  material  strain  hardens  significantly  or 
flows  with  little  hardening,  and  so  on.  The  processes  that  must  be  analyzed  are  inherently 
nonlinear,  but  a  few  general  results  have  been  obtained  for  dynamic  crack  growth  in  elastic- 
plastic  materials.  Some  recent  studies  are  described  in  the  sections  to  follow. 

2.  RATE  INDEPENDENT  MATERIAL  RESPONSE 

Experimental  data  on  the  dependence  of  dynamic  fracture  toughness  versus  crack  speed 
for  A  IS  I  4340  steel  and  other  materials  that  are  commonly  considered  as  elastic-plastic  in 
their  bulk  mechanical  response  have  some  common  features.  Here,  attention  is  limited  to 
situations  in  which  crack  growth  occurs  by  the  single  mechanism  of  void  nucleation  and 
ductile  hole  growth  to  coalescence,  and  in  which  the  extent  of  plasticity  is  sufficiently  limited 


to  permit  interpretation  of  the  fields  surrounding  the  crack  tip  region  on  the  basis  of  a  stress 
intensity  factor.  The  toughness  is  found  to  be  relatively  insensitive  to  variations  in  speed 
for  very  low  speeds  (less  than  20%  of  the  shear  wave  speed)  but  increase  dramatically  with 
crack  speed  for  greater  speeds.  The  speed  dependence  of  the  surrounding  elastic  field  is  not 
nearly  great  enough  to  account  for  this  dependence,  so  an  explanation  must  be  sought  in  the 
plastically  deforming  region  itself.  The  most  likely  reasons  for  this  toughness-speed  behavior 
are  material  inertia  and  material  rate  sensitivity.  While  the  strain  rate  in  the  crack  tip  region 
is  necessarily  very  high,  the  same  general  behavior  has  been  observed  in  materials  that  are 
relatively  rate  insensitive  in  their  bulk  response  up  to  strain  rates  of  10 1  sec-1,  so  the  role 
of  material  inertia  has  been  examined  separately  from  the  role  of  rate  effects  in  research  on 
this  matter.  The  general  idea  is  to  generate  a  theoretical  fracture  toughness  versus  crack 
speed  relationship  to  determine  the  role  of  inertia  on  the  scale  of  crack  tip  plastic  zone  on 
the  observed  dynamic  crack  growth  response. 

A  rough  estimate  of  the  conditions  under  which  material  inertia  has  a  significant  influ¬ 
ence  on  the  development  of  fields  within  the  active  plastic  zone  is  obtained  as  follows.  For 
steady  quasistatic  growth  of  a  crack  in  the  plane  strain  opening  mode  in  an  elastic-ideally 
plastic  material,  Rice,  Drugan  and  Sham  [2]  have  constructed  an  asymptotic  field  consisting 
of  a  constant  state  region  ahead  of  the  crack  tip,  followed  by  a  fan  sector  with  singular  plastic 
strain,  an  elastic  unloading  region,  and  finally  a  small  plastic  re  loading  zone  along  the  crack 
flanks.  Within  the  region  of  singular  plastic  strain,  the  distribution  of  particle'  velocity  and 
shear  strain  in  crack  tip  polar  coordinates  is 

»r,  "0  ~  >'fn  hi  .  fvr8  ~  fn  In  (-^)  (2.1) 

where  f n  is  the  tensile  yield  strain.  r„  is  the  maximum  extent  of  the  plastic  /oik',  and  r  is 


the  radial  distance  from  the  crack  tip.  If  an  expression  for  the  kinetic  energy  density  and 
the  stress  work  density  are  derived  for  this  deformation  field,  then  the  ratio  of  the  kinetic 
energy  density  to  the  stress  work  density  as  a  function  of  r  is 

2 

KE/SW  %  In  (~)  (2.2) 

where  c  =  sj  E  /  p  is  an  elastic  shear  wave  speed.  Based  on  this  simple  estimate,  it  might 
be  expected  that  inertial  effects  will  be  significant  when  the  ratio  in  (2.2)  is  greater  than 
one-tenth.  For  example,  if  v/c  =  0.1  then  the  ratio  is  greater  than  one-tenth  if  r/r„  <  0.3. 

In  retrospect,  this  estimate  has  a  feature  that  could  have  been  interpreted  as  a  warning 
that  asymptotic  analysis  of  this  problem  would  have  some  subtle  difficulties.  The  estimate 
suggests  that  for  any  nonzero  v/c  there  is  a  range  of  r/rQ  for  which  material  inertial  effects 
are  important,  but  that  the  size  of  that  region  diminishes  very  rapidly  as  v/c  approaches 
zero.  Indeed,  the  estimate  suggests  that  inertial  effects  are  important  only  over  a  region  for 
which  (r/r0)('/c  <  1  which  is  similar  in  form  to  the  restriction  on  the  domain  of  validity  of 
the  dynamic  asymptotic  solution  given  below  in  (2.4). 
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The  steady-state  growth  of  a  crack  at  speed  r  in  the  the  antiplane  shear  mode,  or  mode 
III  in  fracture  mechanics  terminology,  under  small  scale  yielding  conditions  was  analyzed  by 
Freund  and  Douglas  [3]  and  by  Dunayevsky  and  Achenbach  [4],  The  field  equations  governing 
this  process  include  the  equation  of  momentum  balance,  the  strain-displacement  relations, 
and  the  condition  that  the  stress  distribution  far  from  the  crack  edge  must  lie  the  same  as 
the  near  tip  stress  distribution  in  a  corresponding  elastic  problem.  For  olastie-ideally  plastic 
response  of  the  material,  the  stress  state  is  assumed  to  lie  on  the  Mises  yield  locus,  a  circle  of 
radius  r„  in  the  plane  of  rectangular  stress  components,  and  the  stress  and  strain  are  related 
through  the  incremental  Prandtl  Reuss  flow  rule.  The  material  is  linearlv  elastic  with  shear 
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modulus  ft  in  regions  when'  the  stress  state  does  not  satisfy  the  yield  condition. 


With  a  view  toward  deriving  a  theoretical  relationsliip  between  the  crack  tip  speed  and 
the  imposed  stress  intensity  factor  required  to  sustain  this  speed  according  to  a  critical  plastic 
strain  crack  growth  criterion,  attention  was  focussed  on  the  strain  distribution  on  the  crack 
line  within  the  active  plastic  zone,  and  on  the  influence  of  material  inertia  on  this  stress 
distribution.  It  was  found  that  the  distribution  of  shear  strain  on  this  line,  say  ■)yz(x,0)  in 
crack  tip  rectangular  .r,  y  coordinates,  could  be  determined  exactly  in  terms  of  the  plastic 
zone  size  r0  in  the  parametric  form 
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where  rn  =  v/c,  and  r_,  is  the  elastic  shear  wave  speed.  While  the  integral  I[t)  has  a 
representation  in  terms  of  elementary  functions  only  for  very  special  values  of  its  argument, 
it  is  easily  evaluated  by  numerical  methods  for  any  nonzero  value  of  m. 

The  exact  result  (2.3)  resolved  a  long  standing  paradox  concerning  mode  III  crack  tip 
fields.  Rice  [2]  showed  that  the  near  tip  distribution  of  strain  7v;(.r,0)  for  steady  growth 
of  a  crack  under  equilibrium  conditions  was  singular  as  liT(.r/r„)  as  x/r„  — *  0.  On  the 
other  hand,  Slepyan  [5]  showed  that  the  asymptotic  distribution  for  any  rn  >  0  was  of  the 
form  (m-1  —  1)  ln(.r/r„)  as  x/rn  — >  0.  These  two  features  could  be  verified  by  examining 
the  behavior  of  the  exact  solution  for  dynamic  growth  (2.3)  under  the  condition  w  — +  0  for 
any  nonzero  value  of  x/r„  and  under  the  condition  that  x/rn  — ►  0  for  any  nonzero  value  of 
rn,  respectively.  The  resolution  of  the  paradox  was  found,  however,  in  the  observation  that 
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Slepyan’s  asymptotic  solution  is  valid  only  if 

,  .  v2m/(l  +  m)  /O  1  \ 

(.r/r0)  <C  1  (2.4) 

Thus,  the  apparent  inconsistency  arises  from  the  fact  that  the  asymptotic  result  clue  to 
Slepyan  is  valid  over  a  region  that  becomes  vanishingly  small  as  m  — ►  0. 

Graplis  of  the  plastic  strain  distribution  on  the  crack  line  in  the  active  plastic  zone  are 
shown  in  Fig.  2.1  for  m  =  0,  0.3.  0.5.  The  plastic  strain  is  singular  in  each  case,  as  has  already 
been  noted.  The  most  significant  observation  concerns  the  influence  of  material  inertia  on 
the  strain  distribution.  An  increase  in  crack  speed  results  in  a  substantial  reduction  of  the 
level  in  plastic  strain  for  a  fixed  fractional  distance  from  the  crack  tip  to  the  elastic-plastic 
boundary.  Therefore,  if  a  local  ductile  crack  growth  criterion  is  imposed,  then  it  would 
appear  that  the  fracture  resistance  or  toughness  would  necessarily  increase  with  increasing 
crack  tip  speed.  To  quantify  this  idea,  the  fracture  criterion  proposed  by  McClintock  and 
Irwin  [6]  was  adopted.  According  to  this  criterion,  a  crack  will  grow  with  a  critical  value 
of  plastic  strain  at  a  point  on  the  crack  line  at  a  characteristic  distance  ahead  of  the  tip. 
The  crack  will  not  grow  for  levels  of  plastic  strain  at  this  point  below  the  critical  level,  and 
levels  of  plastic  strain  greater  than  the  critical  level  are  inaccessible.  To  make  a  connection 
between  the  plastic  strain  in  the  active  plastic  zone  and  the  remote  loading,  a  relationship 
between  the  size  of  the  plastic  zone  and  the  remote  applied  stress  intensity  factor  is  required. 
This  can  be  provided  only  through  a  complete  solution  of  the  problem,  and  it  was  obtained 
for  the  case  of  mode  III  by  Freund  and  Douglas  on  the  basis  of  a  full  field  numerical  solution 
of  the  governing  equations.  The  resulting  theoretical  fracture  toughness  A'///,/  versus  crack 
speed  is  shown  in  Fig.  2.2  for  continuous  variation  of  the  critical  plastic  strain  from  ~r  =  0 
to  .  20 r„/ft.  I  he  critical  distance  has  been  eliminated  in  favor  of  A'///,  ,  the  level  of 
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applied  stress  intensity  required  to  satisfy  the  same  criterion  for  a  stationary  crack  in  the 
same  material  under  equilibrium  conditions.  The  variable  intercept  at  rn  =  0  indicates  an 
increasing  amount  of  plasticity  with  increasing  critical  plastic  strain,  and  the  intercept  values 
correspond  to  the  so-called  steady  state  toughness  values  of  the  theory  of  stable  crack  growth, 
that  is.  with  the  plateau  level  of  the  resistance  curve. 

The  plot  in  Fig.  2.2  illustrates  some  typical  features.  The  ratio  of  Kiud/Kuic  a 
monotonically  increasing  function  of  crack  speed  in  for  fixed  critical  strain,  and  this  function 
takes  on  large  values  for  moderate  values  of  m.  Although  there  is  no  unambiguous  way  to 
associate  a  terminal  velocity  with  these  results,  they  suggest  a  maximum  attainable  velocity 
well  below  the  elastic  wave  speed  of  the  material.  It  is  emphasized  that  the  variation  of 
toughness  with  crack  speed  in  Fig.  2.2  is  due  to  inertial  effects  alone.  The  material  response 
is  independent  of  rate  of  deformation,  and  the  crack  growth  criterion  that  is  enforced  involves 
no  characteristic  time.  If  inertial  effects  were  neglected,  the  calculated  toughness  would  be 
completely  independent  of  speed.  The  question  of  the  influence  of  material  rate  sensitivity 
on  this  relationship  is  a  separate  issue. 

The  equivalent  plane  strain  problem  of  dynamic  crack  growth  in  an  elastic-ideally  plastic 
material  has  not  been  so  fully  developed.  However,  a  numerical  calculation  leading  to  a 
fracture  toughness  versus  crack  speed  relationship,  analogous  to  Fig.  2.2.  has  been  described 
by  Lam  and  Freund  [7].  They  adopted  the  critical  crack  tip  opening  angle  growth  criterion 
and  derived  results  for  mode  I  on  the  basis  of  the  Mises  yield  condition  and  ./•>  flow  theory 
of  plasticity  that  are  quite  similar  in  general  form  to  those  shown  for  mode  III.  The  nature 
of  the  elastic-plastic  fields  deep  within  the  active  plastic  zone  were  difficult  to  discern  from 
the  finite  element  results,  and  an  analytical  study  of  the  asymptotic  field  was  undertaken 
by  Leighton.  Champion  ami  Freund  [S]  in  order  to  examine  this  feature.  It  was  found  that 
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the  plastic  strain  components  had  hounded  limiting  values  at  the  crack  tip  for  any  nonzero 
crack  speed  m  =  e/e,,  hut  that  these  limits  depended  on  crack  speed  as  1/???.  Both  Slepyan 
[5j  and  Achcnhach  and  Dunayevskv  [9]  reported  studies  of  this  problem  in  which  they  took 
elastic  compressibility  into  account.  They  were  able  to  extract  solutions  valid  very  close  to 
the  crack  tip  in  the  limit  of  vanishing  crack  speed.  It  should  be  noted  that  in  both  studies 
the  Tresca  yield  condition  was  used  together  with  the  Mises  flow  rule,  so  that  the  fields 
described  are  consistent  with  normality  of  the  plastic  strain  rate  to  the  yield  surface  only 
in  the  limit  of  incompressibility.  In  addition,  both  studies  imposed  restrictions  on  the  out- 
of-plane  deformation  that  arise  from  assuming  normality  of  the  plastic  strain  rate  to  the 
Tresca  yield  surface.  In  a  study  of  the  same  problem  in  the  incompressible  limit.  Gao  and 
Nemat-Xasscr  [10]  reported  a  solution  with  jumps  in  stress  and  particle  velocity  across  radial 
lines  emanating  from  the  crack  tip  for  all  crack  speeds  between  zero  and  the  elastic  shear 
wave  speed,  where  the  jump  magnitudes  were  subject  to  the  appropriate  jump  conditions. 
It  is  shown  in  [8],  however,  that  if  the  sequence  of  deformation  states  throughout  the  jump 
must  be  admissible  plastic  states,  consistent  with  the  theory  of  mechanical  shocks,  then 
discontinuities  in  the  angular  variation  of  stress  and  particle  velocity  components  around  the 
crack  edge  can  be  ruled  out. 

Some  data  on  the  dynamic  fracture  toughness  of  metals  during  rapid  crack  growth  are 
available.  Rosakis,  Duffy  and  Freund  [11]  used  the  optical  shadow  spot  method  in  reflection 
mode  to  infer  the  prevailing  stress  intensity  factor  during  rapid  crack  growth  in  43-10  steel 
hardened  to  /?<-•  =  45.  This  is  a  relatively  strain  rate  insensitive  material  with  very  little 
strain  hardening,  so  that  the  material  may  presumably  be  modeled  as  elastic- ideally  plastic. 
I  he  observed  toughness  varied  little  with  crack  speed  for  speeds  up  to  about  GOO  to  700  ill /s. 
and  thereafter  the  toughness  increased  sharply  with  increasing  crack  tip  speed.  The  general 
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form  of  the  toughness  versus  speed  data  was  similar  to  the  theoretical  prediction  based  on 
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the  numerical  simulation  reported  by  Lam  and  Freund  [C],  lending  support  to  the  view  that 
material  inertia  on  the  scale  of  the  crack  tip  plastic  zone  has  an  important  influence  on 
the  perceived  dynamic  fracture  toughness.  Similar  data  were  reported  by  Kobayashi  and 
Dally  [12]  who  made  photoelastic  measurements  of  the  crack  tip  stress  field  by  means  of  a 
birefringent  coating  on  the  specimen.  Data  on  crack  propagation  and  arrest  in  steels  were 
reported  by  Dahlbcrg,  Nilsson  and  Brickstad  [13]. 

3.  VISCOPLASTIC  MATERIAL  RESPONSE 

An  estimate  of  the  plastic  strain  rate  near  the  tip  of  an  advancing  crack  may  be  obtained 
as  follows.  Suppose  that  the  yield  stress  in  shear  is  ry  and  that  the  elastic  shear  modulus  is 
//.  so  that  the  yield  strain  is  ry / // .  As  a  rough  estimate  of  the  plastic  strain  rate,  consider 
the  yield  strain  divided  by  the  time  required  for  the  crack  tip  to  traverse  a  region  that  is  the 
size  of  the  active  plastic  zone  at  speed  v.  Following  McClintock  and  Irwin  [G],  if  the  size  of 
the  plastically  deforming  region  is  interpreted  as  the  largest  extent  in  an  elastic  field  of  the 
locus  of  points  on  which  the  maximum  shear  stress  is  ry,  then  the  estimate  of  strain  rate  is 

(iP)est  ^  7vry/^G  (3-l) 

where  energy  release  rate  G  is  the  characterizing  parameter  for  the  elastic  field.  Clearly,  for 
rapid  growth  of  a  crack  in  a  low  toughness  material,  the  strain  rate  estimate  can  be  enormous, 
well  in  excess  of  10*’  sr  r-1 .  Some  results  concerned  with  fairly  large  plastic  strains  (compared 
to  elastic  strains)  under  high  rate'  conditions  are  discussed  in  this  section,  and  the  case  of 
viseoplastic  crack  growth  with  small  plastic  strains  is  considered  in  the  next  section. 

Consider  steady  crack  growth  in  an  elastic-plastic  maternal  for  which  the  flow  stress 
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depends  on  the  rate  of  deformation.  The  particular  material  model  known  as  the  over-stress 
power  law  model  has  been  considered  by  Lo  [14],  Brickstad  [15]  and  a  number  of  other 
authors.  According  to  this  idealization,  the  plastic  strain  rate  in  simple  shear  ~<p  depends  on 
the  corresponding  shear  stress  r  through 

ip  =  it+io{(r  -T,)/fi}n  for  t  >  t,  (3.2) 

where  ~<t  is  the  threshold  strain  rate  for  this  description,  or  the  plastic  strain  rate  when  r  —  t,. 
The  description  also  includes  the  elastic  shear  modulus  //,  the  viscosity  parameter  70.  and 
the  exponent  n.  A  common  special  case  is  based  on  the  assumption  that  the  slow  loading 
response  of  the  material  is  clastic-idcally  plastic  and  that  all  inelastic  strain  is  accumulated 
according  to  (3.2).  For  this  case,  =  0  and  rt  is  the  slow  loading  flow  stress  r„.  For  other 
purposes,  it  is  assumed  that  (3.2)  provides  a  description  of  material  response  only  for  high 
plastic  strain  rates,  in  excess  of  the  transition  plastic  strain  rate  A!t  and  for  stress  in  excess 
of  the  corresponding  transition  stress  level  rf.  For  low  or  moderate  plastic  strain  rates,  the 
variation  of  plastic  strain  rate  with  stress  is  weaker  than  in  (3.2),  and  a  common  form  for 
the  dependence  is  (cf.  Frost  and  Ashby  [1 G] ) 

ip  =  <7i(r)exp{-r72(7-)}  (3.3) 

where  <j\  and  <)2  are  algebraic  functions.  The  marked  difference  between  response  at  low  or 
moderate  plastic  strain  rates  and  at  high  strain  rates  may  be  due  to  a  change  in  fundamental 
mechanism  of  plastic  deformation  with  increasing  rate,  or  it  may  be  a  structure  induced 
transition.  For  present  purposes,  it,  is  sufficient  to  regard  the  difference  as  an  empirical 
observation.  The  two  forms  of  constitutive  laws  (3.2)  and  (3.3)  can  lead  to  quite  different 
results  in  analysis  of  crack  tip  fields  and,  indeed,  the  form  (3.2)  leads  to  fundamentally 


ditr.  ■rent  results  for  ditferent  values  of  1 1  it*  exponent  n. 


Lo  [14]  extended  some  earlier  work  on  the  asymptotic  field  for  steady  quasistatic  crack 
growth  in  an  elastic-viscoplastic  material  by  Hui  and  Riedel  [17]  to  include  inertial  effects. 
In  both  cases,  the  multiaxial  version  of  (3.2)  with  ~,t  —  0  and  r,  -  was  adopted  to  describe 
inelastic  response,  with  no  special  provision  for  unloading.  They  showed  that  for  values  of 
the  exponent  n  less  than  3.  the  asymptotic  stress  field  is  the  elastic  stress  field.  For  values 
of  n  greater  than  3.  on  the  other  hand.  Lo  constructed  an  asymptotic  field  including  inertial 
effects  having  the  same  remarkable  feature  of  complete  autonomy  found  by  Hui  and  Riedel, 
that  is.  it  revealed  no  dependence  on  the  level  of  remote  loading.  For  steady  antiplane  shear 
mode  III  crack  growth.  Lo  found  the  radial  dependence  of  the  inelastic  strain  on  the  crack 
line  ahead  of  the  tip  to  be 

'y. 0)  ~  ("  -  l)(e/y0.r)1/("_,,T,(e/c,)  (3.4) 

where  the  dependence  of  the  amplitude  factor  T /,  on  crack  speed  is  given  grapliically  by 
Lo.  who  also  analyzed  the  corresponding  plant'  strain  problem.  Note  that  as  n  — ►  x  the 
plastic  strain  singularity  becomes  logarithmic.  The  full  field  solution  for  this  problem  under 
small  scale  yielding  conditions  was  determined  numerically  by  Freund  and  Dougins  [IS]. 
The  numerical  results  showed  a  plastic  strain  singularity  much  stronger  than  for  the  rate 
independent  case,  and  it  appeared  from  the  numerical  results  that  the  domain  of  dominance 
of  the  asymptotic  field  within  the  crack  tip  plastic  zone  expanded  with  increasing  crack  tip 
speed.  These  observations  are  consistent  with  (3.4). 
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t.  HIGH  STRAIN  RATE  CRACK  GROWTH 

A  particularly  interesting  class  of  dynamic  fracture  problems  are  those  concerned  with 
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crack  growth  in  materials  that  may  or  may  not  experience  rapid  growth  of  a  sharp  cleavage 
crack,  depending  on  the  conditions  of  temperature,  stress  state  and  rate  of  loading.  These 
materials  may  fracture  by  either  a  brittle  or  ductile  mechanism  on  the  microscale,  and  the 
focus  of  work  in  this  area  is  on  establishing  conditions  for  one  or  the  other  mode  to  dominate. 
The  phenomenon  is  most  commonly  observed  in  ferritic  steels.  Such  materials  show  a  de¬ 
pendence  of  How  stress  on  strain  rate,  and  the  strain  rates  experienced  by  a  material  particle 
in  the  path  of  an  advancing  crack  are  potentially  enormous.  Consequently,  the  mechanics  of 
rapid  growth  of  a  sharp  macroscopic  crack  in  an  elastic- viscoplastic  material  that  exhibits  a 
fairly  strong  variation  of  flow  stress  with  strain  rate  has  been  of  interest  in  recent  years.  The 
general  features  of  the  process  as  experienced  by  a  material  particle  on  or  near  the  fracture 
path  are  straight  forward.  As  the  edge  of  a  growing  crack  approaches,  the  stress  magnitude 
tends  to  increase  there  due  to  the  stress  concentrating  effect  of  the  crack  edge.  The  material 
responds  by  flowing  at  a  rate  related  to  the  stress  level  in  order  to  mitigate  the  influence  of 
the  crack  edge.  It  appears  that  the  essence  of  cleavage  crack  growth  is  the  ability  to  elevate 
the  stress  to  a  critical  level  before  plastic  flow  can  accumulate  to  defeat  the  influence  of  the 
crack  tip.  In  terms  of  the  mechanical  fields  near  the  edge  of  an  advancing  crack,  the  rate  of 
stress  increase  is  determined  by  the  elastic  strain  rate,  while  the  rate  of  crack  tip  blunting 
is  determined  by  the  plastic  strain  rate.  Thus,  an  equivalent  observation  is  that  the  elastic 
strain  rate  near  the  crack  edge  must  dominate  the  plastic  strain  rate  for  sustained  cleavage. 
It  is  implicit  in  this  approach  that  the  material  is  intrinsically  cleavable,  and  the  question 
investigated  in  concerned  with  tire  way  in  which  work  can  be  supplied  to  the  crack  tip  region. 

The  problem  has  been  studied  from  this  point  of  view  by  Freund  and  Hutchinson  [19]. 

They  adopted  the  constitutive  description  (3.2,3)  with  u  =  1.  This  is  indeed  a  situation 
for  which  the  near  tip  elastic  strain  rate  dominates  the  plastic  strain  rate.  Through  an 


approximate  analysis,  conditions  necessary  for  a  crack  to  run  at  high  velocity  in  terms  of 
constitutive  properties  of  the  material,  the  rate  of  crack  growth,  and  the  overall  crack  driving 
force  were  extracted  under  small  yielding  conditions. 

Consider  the  crack  gliding  along  through  the  elastic-viscoplastic  material  under  plane 
strain  conditions.  At  points  far  from  the  crack  edge,  the  material  remains  elastic  and  the 
stress  distribution  is  given  in  terms  of  the  applied  stress  intensity  factor  Kj.  Equivalently, 
the  influence  of  the  applied  loading  may  be  specified  by  the  rate  of  mechanical  energy  flow 
into  the  crack  tip  region  from  remote  points  G ,  and  these  two  measures  are  related  by  means 
of 

G  =  — — — A(v)K]  (4.1) 

where  v  and  E  are  the  elastic  constants  of  an  isotropic  solid  and  .4  is  a  universal  function  of 
the  instantaneous  crack  speed  v.  The  function  has  the  properties  that  .4(0)  =  1,  .I'fO)  =  0 
and  A(v)  — ►  oo  as  v  — ♦  cr.  For  points  near  the  crack  edge  the  potentially  large  stresses  are 
relieved  through  plastic  flow,  and  a  permanently  deformed  but  unloaded  wake  region  is  left 
behind  the  advancing  plastic  zone  along  the  crack  flanks.  For  material  particles  in  the  outer 
portion  of  the  active  plastic  zone  the  rate  of  plastic  straining  is  expected  to  be  in  the  low 
or  moderate  strain  rate  range,  whereas  for  particles  close  to  the  crack  edge,  the  response  is 
modelled  by  the  constitutive  law  (3.2)  with  n  =  1.  Because  of  elastic  rate  dominance,  the 
stress  distribution  within  this  region  has  the  same  spatial  dependence  as  the  remote  field 
but  with  a  stress  intensity  factor  different  from  the  remote  stress  intensity  factor.  The  crack 
tip  stress  intensity  factor,  say  is  assumed  to  control  the  cleavage  growth  process.  The 

influence  of  the  remote  loading  is  screened  from  the  crack  tip  by  the  intervening  plastic  zone, 
and  the  main  purpose  of  the  analysis  is  to  determine  the  relationship  between  the  remote 
loading  and  the  crack  tip  field.  For  present  pmposes,  it  is  assumed  that  the  crack  grows  as 


a  cleavage  crack  with  a  fixed  level  of  local  energy  release  rate,  say  Gcttp.  The  question  then 
concerns  the  conditions  under  which  enough  energy  can  be  supplied  remotely  to  sustain  the 
level  of  energy  release  rate  Gctip  at  the  crack  tip. 

The  matter  of  relating  the  applied  G  to  Grtlp  was  pursued  by  enforcing  an  overall  energy 
rate  balance.  The  balance  may  be  cast  into  the  form 

Gc,,p  =  G  -  -  [  (l A  —  f  U;dy  (4.2) 

r  J A  J-h 

where  ,4  is  the  area  of  the  active  plastic  zone  in  the  plane  of  deformation,  h  is  the  thickness  of 
the  plastic  wake  far  behind  the  crack  tip,  and  U*  is  the  residual  elastic  strain  energy  density 
trapped  in  the  remote  wake.  This  relation  simply  states  that  the  energy  being  released  from 
the  body  at  the  crack  tip  is  the  energy  flowing  into  the  crack  tip  region  reduced  by  the 
energy  dissipated  through  plastic  flow  in  the  plastic  zone,  and  further  reduced  by  the  energy 
trapped  in  the  wake  due  to  incompatible  plastic  strains.  The  expression  is  exact. 

Through  several  approximations,  the  complete  energy  balance  (4.2)  was  reduced  in  [19] 
to  the  simple  form 

G/Gcltp  =  1  +  D(m)Pc  (4.3) 

where  the  dimensionless  parameter  Pc  is  7 oy/ppGctip(l  +  27^/7,,  r<)/3rf3  and  D(m )  is  a 
dimensionless  function  of  crack  tip  speed  m  =  v/cr  and  p  is  the  material  mass  density.  Pc 
is  a  monotonically  increasing  function  of  temperature  for  steels  with  values  in  the  range 
from  about  0  to  10  as  temperature  varies  from  0A’  to  about  400  A".  The  function  D(m)  is 
asymptotically  unbounded  as  m  — *  0  and  m  -*  1,  and  it  has  a  mininnun  at  an  intermediate 
crack  tip  speed.  The  applied  crack  tip  driving  force,  say  Gs,  is  related  to  the  crack  tip  energy 
release  rate  by  G'.y  =  G/(  1  —  v/cr )  for  a  semi-infinite  crack  in  an  otherwise  unbounded  body, 
and  this  relationship  is  adopted  here  as  an  approximation.  A  graph  of  Gs/Grtip  is  shown  in 


Fig.  4.1  in  the  form  of  a  surface  over  the  crack  speed— temperature  plane. 


The  graph  in  Fig.  4.1  gives  the  locus  of  combinations  G.s-,e,T  for  which  steady  state 
propagation  of  a  sharp  crack  can  be  sustained.  The  implication  is  that  if  a  cleavage  crack 
can  be  initiated  for  a  combination  Gs,t\T  that  is  above  the  surface,  then  the  crack  will 
accelerate  to  a  state  on  the  stable  branch  of  the  surface  (i.e.  the  side  with  increasing  Gs 
at  fixed  temperature).  If  the  driving  force  diminishes  as  the  crack  advances,  or  if  the  local 
material  temperature  increases  as  the  crack  advances,  then  the  state  combination  will  move 
toward  the  minimum  point  on  the  surface  at  the  local  temperature.  If  the  driving  force  is 
further  decreased,  or  if  the  temperature  is  further  increased,  then  growth  of  a  sharp  cleavage 
crack  cannot  be  sustained  according  to  the  model.  The  implication  is  that  the  crack  will 
arrest  abruptly  from  a  fairly  large  speed,  and  a  plastic  zone  will  then  grow  from  the  arrested 
crack. 

Of  special  significance  is  the  observation  that,  at  any  given  temperature,  the  variation 
of  required  driving  force  with  crack  speed  has  an  absolute  minimum,  say  G*s/Grtlp.  This 
implies  that,  according  to  this  model,  it  is  impossible  to  sustain  cleavage  crack  growth  at 
that  temperature  with  a  driving  force  below  this  minimum.  Thus,  this  minimum  as  a  function 
of  temperature  may  be  interpreted  as  the  variation  of  the  so-called  arrest  toughness  for  the 
material  with  material  temperature.  This  minimum  is  plotted  against  temperature  for  the 
case  of  mild  steel  in  Fig.  4.2. 

Further  crack  growth  beyond  the  first  arrest  is  possible  if  either  a  ductile  growth  criterion 
can  be  met  or  if  cleavage  can  be  reinitiated  through  strain  hardening  in  the  evolving  plastic 
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zone.  The  details  of  the  model  have  been  refined  through  full  numerical  solution  of  the 
problem  [20],  but  the  essential  features  have  not  changed  with  more  precise  analysis.  A 
modification  of  the  basic  model  was  proposed  by  Mataga,  Hutchinson  and  Freund  [21]  that 


1G 


KL 


fcit J  ^  lv2v  2 if s .a  .vlyl'V-.Ov 


vv 


provides  a  description  that  is  in  bettor  agreement  with  full  field  numerical  solutions  than 
the  model  outlined  above.  In  the  original  development  [20).  it  was  assumed  that  the  plastic 
dissipation  was  completely  controlled  by  the  near  tip  stress  intensity  factor  Held,  say  hi,,,. 
However,  this  stress  intensity  factor,  which  is  asymptotically  correct,  must  give  way  to  the 
far  field  stress  stress  intensity  factor  I\  with  increasing  distance  from  the  crack  edge.  It  was 
observed  in  [21]  that  the  estimate  of  plastic  dissipation  was  improved  significantly,  at  least 
in  comparison  to  finite  element  simulations,  if  the  plastic  dissipation  was  estimated  on  the 
basis  of  a  “mean"  stress  intensity  factor  ^jI\tipK .  The  graph  of  arrest  toughness  versus 
temperature  for  parameters  corresponding  to  mild  steel  are  also  shown  in  Fig.  4.2. 


Important  experiments  on  crack  propagation  and  arrest  in  steel  specimens  are  currently 
being  carried  out  by  deWit  and  Fields  [22].  Their  specimens  are  enormous  single  edge  notched 
plates  loaded  in  tension.  The  growing  crack  thus  experiences  an  increasing  driving  force  as 
it  advances  through  the  plate.  A  temperature  gradient  is  also  established  in  the  specimen 
so  that  the  crack  grows  from  the  cold  side  of  the  specimen  toward  the  warm  side.  Based 
on  the  presumption  that  the  material  becomes  tougher  as  the  temperature  is  increased,  the 
crack  also  experiences  increasing  resistance  as  it  advances  through  the  plate.  The  specimen 
material  is  A533B  pressure  vessel  steel,  which  is  both  very  ductile  and  strain  rate  sensitive. 
In  the  experiments,  the  fracture  initiates  as  a  cleavage  fracture  and  propagates  at  high  speed 
through  the  specimen  into  material  of  increasing  toughness.  The  crack  then  arrests  abruptly 
in  material  whose  temperature  is  above  the  nil  ductility  temperature  for  the  material  based 
on  (harpy  tests.  A  large  plastic  zone  grows  from  the  arrested  crack  edge,  and  cleavage 
crack  growth  is  occasionally  reinitiated.  The  essential  features  of  the  experiment  appear 
to  be  consistent  with  the  model  of  high  strain  rate  crack  growth  outlined  in  [19],  .and  this 


model  appears  to  provide  a  conceptual  framework  for  interpretation  of  the  phenomenon.  An 
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analysis  of  rapid  crack  growth  in  a  rate  dependent  plastic  solid  lias  also  been  carried  out  by 
Brickstad  [15]  in  order  to  interpret  some  experiments  on  rapid  crack  grou  t  h  in  a  high  strength 
steel.  It  is  note'll  that  the  work  of  deWit  and  Fields  [22]  is  part  of  the  Nuclear  Regulatory 
Commission's  Heavy  Section  Steel  Technology  Program  which  sustains  an  integrated  dynamic 
fracture  mechanics  effort  involving  experiments,  computation  and  mate-rial  characterization. 


5.  CONCLUDING  REMARKS 

The  results  described  in  the  preceding  sections  reflect  some  progress  toward  discover)-  of 
the  role  played  played  by  crack  tip  plastic  fields  in  establishing  conditions  for  rapid  advance 
of  a  crack  in  an  elastic-plastic  material.  Understanding  of  this  issue  is  far  from  complete,  and 
a  few  of  thi'  open  questions  that  could  be  profitably  pursued  are  identified  in  this  concluding 
section.  For  example,  much  of  the  modeling  that  has  resulted  in  a  detailed  description  of 
crack  tip  elastic-plastic  fields  has  been  based  on  the  assumption  that  the  fields  are  steady  as 
seen  by  a  crack  tip  observer.  This  approach  overlooks  all  transient  aspects  of  the  process. 
The  picture  of  the  way  in  which  a  crack  tip  plastic  zone  develops  in  a  cracked,  rate  sensitive 
structural  material  under  the  action  of  stress  wave  loading  is  not  clear,  but  the  question  is 
important  in  the  sense  that  these  fields  determine  whether  or  not  the  crack  will  advance. 
The  same  issue  appears  to  be  at  the  heart  of  the  cleavage  initiation  process  in  steels,  but  on 
a  microstructural  scale.  Here,  the  sudden  cracking  of  carbides  or  other  brittle  phases  due  to 
incompatible  plastic  strains  provides  a  nucleation  mechanism,  and  the  question  is  whether 
or  not  these  dynamic  microcracks  penetrate  into  the  adjacent  ferrite  as  sharp  cracks.  The 
answer  seems  to  hinge  on  the  way  in  which  plastic  strains  develop  near  the  carbide-ferrite 
interface  due  to  the  appearance  of  the  microcracks  in  the  brittle  phases. 

The  transients  of  the  arrest  of  a  cleavage  crack  in  a  structural  material  are  also  unclear 
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at  this  point.  A  running  crack  appears  to  arrest  because  conditions  for  the  continuous 
reinitiation  of  cleavage  cannot  be  maintained  [-31.  In  terms  of  the  model  discussed  in  section 
4.  arrest  occurs  because  conditions  for  elastic  rate  dominance  of  t lit-  local  field  cannot  be 
maintained.  However,  the  model  does  not  provide  information  on  the  process  thereafter.  It 
appears  from  the  experiments  reported  by  deWit  and  Fields  [22]  that  arrest  is  quite  abrupt, 
that  a  largo  plastic  /one  grows  from  the  crack  edge  following  arrest  of  the  cleavage  crack,  that 
the  crack  may  grow  subsequently  in  a  ductile  mode,  and  that  cleavage  may  be  reinitiated 
at  a  later  stage.  It  is  not  clear  if  the  cleavage  reinitiation  is  due  to  a  rate  effect  or  to  a 
combination  of  strain  hardening  and  constraint  in  the  interior  portions  of  tin*  specimen. 

Modeling  of  plasticity  effects  in  dynamic  crack  growth  has  been  restricted  to  two  dimen¬ 
sional  systems,  for  the  most  part.  It  is  likely  that  a  number  :>f  three  dimensional  effects  arc  of 
sufficient  importance  to  warrant  further  investigation.  For  example,  crack  propagation  stud¬ 
ies  are  often  carried  out  with  plate  specimens.  For  such  specimens,  t lit'  transition  from  plane 
stress  conditions  in  regions  far  from  the  crack  tip  compared  to  plate  thickness  to  plane  strain 
or  generalized  plant'  strain  conditions  near  the  crack  edge  is  not  clear.  \ang  and  Freund  [24] 
suggest  that  plane  stress  conditions  prevail  only  for  points  beyond  about  one-half  the  plate 
thickness  from  the  crack  edge  for  elastic  deformations.  Out-of-plane  inertia  is  of  potential 
importance  in  these  three  dimensional  fields,  but  this  effect  has  not  been  investigated  to 
date.  Furthermore,  the  role  of  ductile  shear  lips  at  the  free  surfaces  or  of  ductile  ligaments 
left  behind  a  cleavage  crack  as  it.  advances  through  a  structural  metal  are  not  dear  at  this 
time.  In  a  study  of  fracture  initiation  in  dynamically  loaded  specimens  of  a  ductile  material 
by  Nakamura.  Shill  and  Freund  [25],  it  was  shown  that  these  three  dimensional  effects  are 
potentially  very  significant. 

\  irfuallv  all  of  the  foregoing  discussion  has  focussed  on  issues  from  the  traditional 


framin'  mechanics  point  of  view.  that  is.  the  use  of  a  single  parameter  to  characterize  the 
mechanical  state  of  a  dominant  crack  in  a  stressed  body.  This  section  is  concluded,  however, 
with  mention  of  somewhat  speculative  interpretations  of  recent  work  that  suggest  a  departure 
from  the  traditional  fracture  mechanics  viewpoint.  Some  exciting  new  data  on  fracture 
initiation  and  crack  growth  in  a  4340  steel  in  a  very  hard  condition  were  recently  reported 
by  Ravichandran  and  Clifton  [26].  Through  a  modification  of  the  plate  impact  apparatus, 
they  were  able  to  examine  fracture  initiation  and  crack  growth  of  a  few  millimeters  for  the 
plane  strain  situation  of  a  semi-infinite  crack  in  an  unbounded  body  subject  to  plane  wave 
loading,  at  least  for  a  microsecond  or  two.  At  the  lowest  testing  temperature  reported,  the 
cracks  grew  as  cleavage  cracks.  Based  on  optical  measurements  of  the  surface  motion  of  the 
specimen  and  comparison  with  detailed  elastic-viscoplastic  calculations,  it  appeared  that  the 
cracks  grew  more  nearly  at  constant  velocity  crack  than  with  a  fixed  level  of  energy  releast' 
rate  or  stress  intensity  factor.  This  observation  is  similar  to  that  made  by  Ravi-Chandar  and 
Knauss  [27]  who  studied  crack  growth  in  the  brittle  polymer  Homalite- 100.  In  both  cases, 
this  observation  was  made  in  situations  where  the  load  was  suddenly  applied  and  the  load 
level  was  very  intense  compared  to  the  minimum  load  necessary  to  induce  fracture  in  the 
same  situation.  The  results  suggest  that  the  one  parameter  characterization  of  the  crack  tip 
conditions  may  not  be  adequate  to  describe  fracture  response  under  such  conditions. 
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FIGURE  CAPTIONS 


Fin.  2.1  Total  strain  on  the  crack  line  in  the  active  plastic  zone  for  steady  dynamic  growth 
of  a  mode  III  crack  in  an  elastic-ideally  plastic  material  from  (2.3). 

Fig.  2.2  Theoretical  fracture  toughness  versus  crack  speed  for  steady  growth  of  a  mode  III 
crack  according  to  the  critical  plastic  strain  at  a  characteristic  distance  criterion, 
for  levels  of  critical  plastic  strain  from  7<-  =  0  to  7<-  =  20 r„///..  The  plateau  is  an 
artifice  arising  from  truncation  of  the  surface  at  a  suitable  level. 

Fig.  4.1  A  surface  representing  conditions  on  applied  crack  tip  driving  force  G's,  crack  tip 
speed  r  and  temperature  that  correspond  to  s  The  plateau  is  an  artifice  arising  from 
truncation  of  the  surface  at  a  suitable  level,  toady  crack  propagation,  as  predicted 
by  (4.3). 

Fig.  4.2  The  minimum  driving  force  G'*-  needed  to  drive  the  crack  dynamically  as  a  function 
of  temperature.  The  solid  line  corresponds  to  Fig.  4.1.  and  the  dashed  line  from  a 
modified  result  in  [21]. 
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Fig.  2.1  Total  strain  on  the  crack  line  in  the  active  plastic  zone  for  steady  dynamic  growth 
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for  levels  of  critical  plastic  strain  from  yc  =  0  to  yc  =  20 r0//i.  The  plate; 


'an  is  an 


artifice  arising  from  truncation  of  the  surface  at  a  suitable  level. 

Fig.  4.1  A  surface  representing  conditions  on  applied  crack  tip  driving  force  Gs-  crack  tip 
speed  r  and  temperature  that  correspond  to  steady  crack  propagation,  as  predicted 
by  (4.3).  The  plateau  is  an  artifice  arising  from  truncation  of  the  surface  at  a 
suitable  level. 

Fig.  4.2  The  minimum  driving  force  6’*-  needed  to  drive  the  crack  dynamically  as  a  function 
of  temperature.  The  solid  line  corresponds  to  Fig.  4.1,  and  the  dashed  line  from  a 
modified  result  iu  (21|. 
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